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The relation between relaxation and diffusion is investigated in a Hamiltonian system of globally 
coupled rotators. Diffusion is anomalous if and only if the system is going towards equilibrium. 
The anomaly in diffusion is not anomalous diffusion taking a power-type function, but is a transient 
anomaly due to non-stationarity. Contrary to previous claims, in quasi-stationary states, diffusion 
can be explained by a stretched exponential correlation function, whose stretching exponent is almost 
constant and correlation time is linear as functions of degrees of freedom. The full time evolution is 
characterized by varying stretching exponent and correlation time. 
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I. INTRODUCTION 

Relaxation to thermal equilibrium has been studied in 
Hamiltonian systems with long-range interactions 0, 0, 
H) Q IE IS ' One of the characteristic phenomena in the 
relaxation process is anomalous diffusion, since normal 
diffusion is expected at equilibrium. Anomalous diffusion 
was firstly investigated in a one-dimensional chaotic rnap 
to describe enhanced diffusion in Josephson junctions H, 
and is observed in many systems both numerically 0,13 
UmilJ and experimentally [l^. 

Anomalous diffusion is also observed in Hamiltonian 
dynamical systems. It is explained as due to power-type 
distribution functions U ITsL IT3 | of trapping and untrap- 
ping times of the orbit in the self-similar hierarchy of 
cylindrical cantori Self-similarity is expected to be 
one of the important concepts to understand statistics 
and motion in Hamiltonian systems, but cannot be the 
main feature in systems with many degrees of freedom. 
Then, as the first step of approaching the study of self- 
similarity, we have to clarify when anomalous diffusion 
appears, and what is the origin of the anomaly. 

Latora et al. 0| discussed the relation between the 
process of relaxation to equilibrium and anomalous dif- 
fusion in a globally coupled rotator system, by compar- 
ing the time series of the temperature and of the mean 
squared displacement of the phases of the rotators. They 
showed that anomalous diffusion changes to a normal dif- 
fusion after a crossover time, and that the crossover time 
coincides with the time when the canonical temperature 
is reached. They also claim that anomalous diffusion oc- 
curs in the quasi-stationary states, which appear before 
the system goes towards equilibrium. 

The crossover from anomalous to normal diffusion de- 
termines the time when the anomalous diffusion finishes. 
However, it is not clearly pointed out when the anoma- 
lous diffusion starts, and hence the study of the relation 
between the relaxation process and anomalous diffusion 



is still not complete. Moreover, in Ref.[lOj, the numeri- 
cal calculations were performed by using only one type of 
initial condition, but different types of initial condition 
may change the conclusion . 

In this article, we study the globally coupled rotator 
system considered in Ref. Iff] , and we exhibit the relation 
between relaxation to equilibrium and anomalous diffu- 
sion with a different type of initial condition from one 
used in Ref.^3|- Then we show that diffusion becomes 
anomalous if and only if the state is neither stationary nor 
quasi-stationary. In other words, diffusion is shown to be 
normal in quasi-stationary states, although a stretched 
exponential correlation function is present, contrar y to 
previous claims that report power law type function |l7l| . 
Simple scaling laws of the correlation function imply that 
the result holds irrespective of degrees of freedom. 

This paper is organized as follows. The model, initial 
condition and observed quantities are described in SeclTTl 
In Sec lIIII we study relaxation process, which we divide 
into three stages, quasi-stationary, relaxational and equi- 
librium stages. Diffusion process in each stage is inves- 
tigated in Sec lIVI by using stretched exponential corre- 
lation functions of momenta. Dependence on degrees of 
freedom is also reported both in Secs lIIII and IIVI The 
last sectionals devoted to summary. 



II. MODEL, INITIAL CONDITION AND 
OBSERVED QUANTITIES 

The model considered in this paper has N classical and 
identical rotators confined to move on the unit circle, and 
the Hamiltonian is composed of a kinetic and a potential 

part (aaiiiaiia , 
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The N particles are globally coupled through the mean 
field defined as 

1 ^ 

M = — ^(cos6lj,sin6'j) = M(cos 0, sin <?!)) , (2) 

3 = 1 

where the modulus M (0 < M < 1) represents the mag- 
netization of this system. We remark that the potential 
V and the kinetic energy K are related to the magneti- 
zation M as follows: 



2V/N = 1- M^, 2K/N = 2U-1 + M^ 



(3) 



where U is the energy per particle, i.e. U = E/N, and 
E is the total energy. The free energy of this system has 
been obtained in the canonical ensemble Jj, , and 

it has been shown that system has a second-order 
phase transition at the critical energy Uc = 0.75. If the 
energy U is greater than the critical energy, the largest 
Lyapunov exponent goes to zero in the thermodynamic 
limit (N oo) Then, all rotators freely rotate, 

and diffusion becomes ballistic. On the contrary, if U 
is small compared to Uc, all rotators are trapped in the 
potential well and no diffusion occurs. We are therefore 
interested in a value of the energy which is near but less 
than the critical energy in order to allow some particle 
diffusion. Hereafter, we set U = 0.69 (a value studied 
also in Refs.[ll[l3). 

The canonical equations of motion for system can 
be cast in a form that uses the mean field Q as follows: 



time series of the temperature, since the canonical aver- 
age of 2K/N coincides with the canonical temperature. 

To observe the diffusion process, we introduce the 
mean square displacement of phases (t) defined as 



1 ^ 

i=i 

The symbol represents the average over all the N 
rotators. The quantity cr|(t) typically scales as ag{t) ~ 
t", and the diffusion is anomalous when a ^ 1,2, while 
it is normal when a = 1 and ballistic for a = 2. The 
quantity crg(t) can be rewritten by using the correlation 
function of momenta Cp(t; r) as 



Ul(t) = / dtl / At2 {Pj{tl)Pjit2))M 

Jo Jo 

pt nt — S 

= 2 ds dTCp{s;T), 
Jo Jo 



(6) 



where Cp{t;T) is defined as 



N 



(7) 



Moreover, if the system is stationary and Cp{t;T) does 
not depend on r accordingly. 



Cpit;T) = Cp{t;0), (Vr>0) 
then Eq.© is simplified as 



(8) 



d^ 

dt 



Pv ^ = -M{t) sm{9,-<P{t)), (j = 1, • • • , N). 

(4) 

We numerically integrate Eq. by using 4th order sym- 
plectic integrators 20^ ^] . The time slice of the integra- 
tor is set at At = 0.2 or 0.4, and it suppresses the relative 
energy error down to \AE/E\ < 5 x 10~^. 

We have performed the integrations starting from 
M{0) — 0. To prepare these initial conditions numeri- 
cally, we set (?j(0) = 2'Kj/N , and Pj{G) is taken from a 
uniformly random distribution whose support is [—p,p\, 
where the value p is chosen to get the energy density 
U. The total momentum ^J^iPj is an integral of the 
motion and we initially set it to zero. This initial state 
corresponds to a local entropy minimum 22], and to a 
stationary stable solution to the Vlasov-Poisson equation 
, although the systenigoes towards Gibbs equilibrium 
due to finite size effects 'q. With respect to the M(0) = 1 
initial condition chosen in Ref . (IQj , the one we choose has 
the advantage of being a quasi-stationary state from the 
start. 

We numerically observe the time series of two quanti- 
ties. One is for the relaxation process, and the other is 
for the diffusion process. 

To observe the relaxation process, we use the magne- 
tization M{t). Note that observing M{t) corresponds to 
observing 2K{t)/N by using Eq.Q, and 2K{t)/N is the 



al{t) = 2 / {t-s) Cp{s-0) ds. (9) 
Jo 



III. RELAXATION PROCESS 

Temporal evolutions of M{t) are shown in FigQ] In 
order to suppress fiuctuations, we have calculated aver- 
ages over realizations. Throughout this article, unless 
no comments appear, the numbers of realizations are 
n = 1000, 100 and 8 for iV = 100, 1000 and 10000, re- 
spectively. We divide the temporal evolutions into three 
stages, Stage-1,11 and III. In Stage-I, the value of magne- 
tization is almost constant but smaller than the canonical 
value. After Stage-I, magnetization rapidly increases to- 
wards its equilibrium value Moq, and we call this time 
interval Stage-II. Finally the system reaches equilibrium, 
during Stage-Ill. 

Let us define boundary times between Stage-I and II, 
ti/iij a-nd between Stage-II and III, tn/m, as follows. 
The magnetization takes the local minimum at i niin i and 
we adopt ij/n = ^min- We define the other bound- 
ary time tii/iii as the first passage time which satisfies 
M{t) = 0.99Moq. Values of the two boundary times are 
reported in Fig|2las functions of degrees of freedom. The 
local minimum time is proportional to A^^'^ for N > 100 
with our initial condition M(0) = 0, as with another 
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FIG. 1: Temporal evolutions of M{t). U = 0.69 and N = 
100, 1000, 10000. The horizontal line represents the canonical 
equilibrium value of M. On each curve, two short vertical 
lines are marked. The first and the second ones are at the 
end of Stage-I and II, respectively. Solid curves are hyperbolic 
tangent functions llOI I. 



initial condition M(0) = 1 |22.j. For small N, we can- 
not neglect the initial time region i < 6 in which the 
level of M{t) goes to 0{1/^/N) coming from the law of 
large numbers (see Fig|2{d)), and hence the power law 
breaks. The power law recovers by subtracting the ini- 
tial increasing time 6 from ti/n as shown in FigEl i-C. 



6 - N^-^ (N > 10). 




FIG. 2: Dependence on degrees of freedom of tj/u (squares) 
and tii/iii (crosses). Stars represent ti/n — 6. The lower 
straight line represents the power law Ai'^'^/150. The upper 
curve is a theoretical prediction of the boundary time tn/m 
using Eqs.JTTJ and JHJ with Mth = 0.99Afcq. 



ac = dM/d(logjo *)|iogio t=bi and lO'' is the time scale. 
As shown in FigO these four parameters are fitted as 



a{N) = 



lO^W = iiVl■^ 

lOOc(iV)' 9 



^(^)^(Me,-d(iV))^ ^^^^ 



1.7 



(11) 



By using the scaling law (|ll|l . we can predict when M{t) 
reaches a given threshold level, A/th, as a function of N. 
Let tth be the threshold time, which satisfies M{tth) = 
Mth, then tth is expressed as 



tth = 10' 



th 



ilfeq - Mth 



In 10 

2a 



(12) 



In Fig|21 tth is reported for Mth — 0.99Meq, and the 
prediction is in good agreement with numerical results. 
We remark that, roughly speaking, ti/n is asymptotically 
proportional to N-^-'^ . 






FIG. 3: Four parameters a, b, c and d are reported as func- 
tions of degrees of freedom, (a) Log-log plot of ac. (b) Log-log 
plot of 10*". (c) Liner- log plot of c. (d) Log-log plot of d. Solid 
curves are scaling functions described in Eg. llH . 



The system seems quasi-stationary in Stage-I. The ex- 
istence of quasi-stationary states for sufficiently long time 
has been questioned in Ref.[23. We will answer to the 
question by observing dependence on r of the correlation 
function Cp{t;T) in Sec lIV Bl 



A theoretical prediction of in/m, the upper curve in 
Fig|21 is obtained by fitting the magnetization M{t) as 
hyperbolic tangent function, 

M {t) = [1 + tanh(a(logio t-b))] c + d. (10) 

The parameter d represents the initial level of M{t), and 
c the half width between initial and equilibrium levels 
of M{t). The product ac is the slope at logj^g t = 6, i.e. 



IV. DIFFUSION PROCESS 

As described in Eq.©, the mean square displacement 
CTg(t) is obtained from correlation function of momenta 
Cp{t\ r), and hence we study diffusion process by observ- 
ing the correlation function. We start from the simplest 
stage, Stage-Ill, because we may use the simple expres- 
sion l^. Next, we progress to Stage-I, where we expect 
that the system is quasi-stationary, and that we may use 
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Eq.(|5J again. In non-stationary stage, Stage-II, we check 
whether diffusion is of a power type. FinaUy we investi- 
gate dependence on degrees of freedom for some impor- 
tant parameters. 



A. Diffusion at Equilibrium 

Assuming the system has reached equihbrium at t 
toq, we observe Cp{t;teq) and ag{t;tcq), where 



<^Ut;toq)^{[0,{t + t,q)~9,{t,q)f) 



N 



(13) 



At equihbrium we may assume that the system is sta- 
tionary, 



Cpit;Uq + T) = Cp{t;t,q), (Vt>0) (14) 



and hence 



alit; icq) = 2 / (t-s) Cp{s; t,q) ds. (15) 
Jo 

Now let us consider the correlation function for N = 
1000. We adopt Uq = 2^° ~ 10^ which is long enough 
to reach equihbrium imaging from Fig^ The correlation 
function Cp(i;teq) is reported in Fig^fa), and is well ap- 
proximated by the stretched exponential function p3 | 



Cp{t;t,q) = 0.47exp[-(i/410)'^-^1 



(16) 



rather than by a pure exponential (see the inset of 
Fig^a) which is a log-linear plot of Cp{t; icq))- 

We remark that a stretched exponential function 
exp[— x'^] with a small exponent |/3|<<cl is indistinguish- 
able from a power- type function in the region |/31na;|-^l: 

exp[— x^] = exp[— exp(/31n2;)] '--^ exp[— 1— /31nx] = x^^ /e. 

However the fitting function Hlt)|) well agree with numeri- 
cal result even around |0.32 ln(t/410)| = 1, whose two so- 
lutions are t ~ 18, 9330. We therefore adopt a stretched 
exponential function as an approximation of Cp{t] tcq). 

By using the fitting function (|16|l and Ea. H15|l . we 
numerically reproduce a'^{t;teq), and the reproduced 
curve well approximates the numerical result as shown 
in Fig0{b). Note that (Tg(i;teq) is proportional to in 
the limit of i ^ 0, since Cp{s\tcq) in Ea. (|15|) goes to 
the constant Cp (0; icq)- On the other hand, in the limit 
of t ^ oo, a1(t;toq) is proportional to t, because both 
Cp(s;teq) and sCp{s;tcq) are almost zeros in long time 
region, and hence their integrals become constants. The 
crossover from to t is also observed if we assume an 
exponential correlation function, and hence we conclude 
that diffusion at equilibrium is normal as expected al- 
though a stretched exponential is present. 



Fluctuation Level 
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FIG. 4; (a) Double log-log plot of normalized correlation 
function Cp(t; feq)/Cp(0; teq) at equilibrium with icq ~ 2^". 
A'^ = 1000. We take an average over n = 100 realizations. The 
straight line and the curve represent the stretched exponen- 
tial function llbll and the power- type function (f/410)"°'^Ve, 
respectively. The upper horizontal line is fluctuation level 
0{l/VNn). The inset shows a log-linear plot with the 
stretched exponential function, (b) Log-log plot of ag{t;tEq) 
with the approximate function produced by Egs. 11511 and 1161 . 



B. Diffusion in Quasi-Stationary State 

Except for Stage-Ill, we cannot expect stationarity to 
hold (jS)) any more. However, from the temporal evolu- 
tions of M{t), Fig^ we may expect quasi-stationarity in 
Stage-I, 



Cp{t-T)=Cp{t-Q) + e{t-T), 



(17) 



where r belongs to Stage-I and e{t) is suitably small. 

The correlation function Cp{t]T) for various values of 
T is reported in Figsl^a) and (b) for N — 1000, and the 
relative error of correlation function, defined as 



Rir) 



mT)\ 

max — 
t Cp(0;0)' 



(18) 



is also reported in Fig[SJc) as a function of t. The error 
i?(r) stays small up to the end of Stage-I, and hence we 
conclude that the system is quasi-stationary in Stage-I. 
We believe that the quasi-stationary states correspond to 
stationary stable states of the Vlasov equation ||25| . We 
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FIG. 5: Correlation function of momenta Cp{t; r) for various 
values of r = 0, 1024, 2048, 4096, 16384, 65536 and 1048576 
from left to right, (a) Log-log plot, (b) Double log-log plot, 
(c) Relative error llTSt of Cp(t;r) from Cp(t;0). 
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FIG. 6: Correlation function of momenta at r = 0, i.e. 
Cp{t;0). The inset is magnification of the horizontal axis 
around t = for N = 100. These numerical results are ap- 
proximated by solid curves which are stretched exponential 
functions 11911 . 
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FIG. 7: Time series of the mean square displacement of the 
phases (jl{t). N = 100, 1000 and 10000 from top to bottom. 
The vertical axis is the original scale only for A'' = 10000, and 
is muhiplied by 10^ and lO*' for iV = 1000 and 100 respectively 
just for a graphical reason. In Stage-I where the system is 
quasi-stationary, the numerical results are approximated by 
solid curves which are obtained from Eq.@ using functions 
H19|l . After the system reaches equilibrium, difi^usion becomes 
normal. Anomaly in diffusion is observed only in Stage-II. 
The two short vertical lines on each curve show the end of 
Stages-I and II, which correspond to the ones found in FigQ 



remark that R{t) is constant in Stage-Ill again due to 
stationarity at equilibrium. 

It seems natural that we regard Cp(t; t) as a series of 
stretched exponential functions of t rather than power- 
type functions, since this function fits Cp{t]T) in more 
than two decades of time (power law fits of the correlation 
functions hold in one decade). Moreover, at equilibrium, 
Cp(t;tcq) is also a stretched exponential rather than a 
pure exponential, as shown in Fig0| 

In the quasi-stationary region, Stage-I, the mean 
square displacement Cg^t) can be derived by the corre- 
lation function Cp(t;0), which is reported in FigEl for 



N = 100, 1000 and 10000. We approximate Cp(i; 0) by a 
stretched exponential function as 

TV = 100 : Cp{t; 0) = 0.38 exp[-(t/20)" *^^], 
N = 1000 : Cp{t; 0) = 0.38 cxp[-(t/180)°-^^], (19) 
N = 10000 : Cp{t; 0) = 0.38 exp[-(t/2200)°-^°]. 

The prefactor 0.38 comes from Cp(0; 0) = 2K{0)/N. 

Using the approximate functions ((T^ and Eq.®, we 
are able to reproduce crg{t), as shown in Figd The 
approximation is good in Stage-I, i.e. in the quasi- 
stationary time region, irrespective of the value of N. 
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Consequently, there is no anomaly in diffusion in Stage- 
I, since the diffusion is explained by stretched exponential 
correlation function. 

C. Diffusion in Non-Stationary State 

After the quasi-stationary region, diffusion becomes 
anomalous, which is faster than normal diffusion, in 
Stage-II. If we fit <Tg{t) by a power- type function t" in 
Stage-II, the exponent a is estimated as 1.54, 1.59 and 
1.74 for N = 100, 1000 and 10000 respectively. The 
values of exponent tend to increase as TV increases as re- 
ported for the systein_having the so-called 2-dimensional 
egg-crate potential 3]. On the other hand, the dura- 
tion in which diffusion is anomalous becomes shorter and 
shorter in logarithmic time scale as N increases , in ac- 
cordance with the sharper change of M(t). Moreover, 
(Tg/f^ is not constant, but has a wave in Stage-II (see 
Fig|SJ. Hence we guess that the anomaly in diffusion 
is not anomalous diffusion taking a power-type function, 
but a transient anomaly due to non-stationarity of Stage- 
II. 
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FIG. 8: Log- log plot of ag{t)/t°'. The exponent a is esti- 
mated as 1.54, 1.59 and 1.74 for = 100, 1000 and 10000 
respectively. The two short vertical lines on each curve show 
the end of Stages-I and II. In Stage-II, o-g/t" is not constant. 
The vertical axis is the original scale only for A'^ = 100, and 
is multipled by 10 and 100 for N = 1000 and iV = 10000 
respectively for a graphical reason. 

Let us proceed to investigate the origin of anomaly in 
diffusion. We focus on the behavior for N — 1000. The 
mean square displacement ag{t) is perfectly determined 
by the correlation function Cp{t;T) using Eq.®, once 
we assume that Cp{t; r) is a series of stretched exponen- 
tial functions. We introduce three parameters, Cp(0;T), 
tcovriT) and /3(t), to describe the stretched exponential 
function as 

Cp{t;T) = Cp(0;T)exph{VUrr(T)}'^(")]. (20) 

We investigate which of the three parameters is the most 
important to yield anomaly in diffusion. 



The strategy is as follows. We reproduce d(jg{t)/dt by 
using the three parameters and the formula 

^(i) = 2 J drCp(0;r)exp[-{(t-r)/W.(r)}'^M]. 

(21) 

We consider the first derivative of ag instead of ag it- 
self, because the former requires only single integration 
while the latter requires double integrations 10. We first 
omit the dependence on r of the parameter Cp(0; t) and 
fix it to a constant value to observe how it affects the 
anomaly in diffusion. We then fix the two other param- 
eters tcoiriT) and /3(r) to determine their effect on the 
mean square displacement. 

From the numerical results of Cp{t; r), Fig|SJb), we de- 
termine the values of three parameters Cp(0; r), tcoiT(''') 
and /3(r) at some value of r by using the least square 
method. The discrete values of the parameters are not 
enough to reproduce dag{t)/dt accurately, and then we 
approximate the parameters by hyperbolic tangent func- 
tions as follows: 

Cp(0; t) = 0.046 [1 + tanh(2.5(logio r - 4.35))] + 0.385, 
tcorr(r) 80 [1 + tanh(1.5(logio t - 3.4))] + 170, 
/3(r) = 0.31 [I + tanh(1.5(logio r - 3.8))] + 0.29. 

(22) 
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FIG. 9: The three parameters Cp(0;r), tcorr(T) and /3(r) as 
functions of r. The latter two parameters tcorr(T) and /3(r) are 
multiplied by 1/1000 and 1/2 respectively for a graphical rea- 
son. Solid curves are hyperbolic tangent functions described 
in Ea. lf^ . 

The hyperbolic tangent functions are in good agree- 
ment with numerical results, as shown in FiglHl To con- 
firm the validity of the approximation, we reproduced 
dag/dt using Eas. (|2I() and (|22|l . and the reproduced one 
is in good agreement with numerical results, as shown in 
FigCnta). 

If we fix Cp(0; t) at its middle value, 0.431, we find that 
the dependence on r of Cp(0;T) does not affect signifi- 
cantly dag/dt, as shown in Fig llOf b). By fixing icorr(''') 
at its middle value, 250, we obtain the same conclusion 
for tcorr(T) as for Cp(0;r), particularly in Stage-II (see 
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Fig llOf c')'). On the contrary, if we fix (3{t) at 0.6 or 0.9, 
we observe no anomaly in diffusion as shown in Fig llOf d). 
because dag (t) /dt is proportional to t and is constant in 
short and long time regions respectively, and the same 
behavior is obtained at equilibrium (see Fig^fb)). Con- 
sequently, among the three parameters, /3(t) plays a cru- 
cial role to produce anomaly in diffusion. 






FIG. 10; Time derivative of the mean square displacement , 
dag{t)/dt. (a) Numerical results (crosses) and reproduced one 
(solid curve) using Ea. l2H and the approximate functions of 
the three parameters in Eas. l22ll . In (b), (c) and (d), Cp(0; r), 
icorr(T") and /3(t) are kept constant, respectively. In (d), two 
constants for /3(t) have been tested. Solid and dashed curves 
represent f3 = 0.6 and 0.9, respectively. The short vertical 
hnes mark the end of Stages-I and II. 



D. Dependence on Degrees of Freedom 




FIG. 11: Double log-log plots of correlation functions for var- 
ious values of degrees of freedom A'^. (a) Cp{t;0) (Stage-I). 
(b) Cp{t-Uq) (Stage-Ill). In both (a) and (b), iV = 100(1000), 
200(500), 300(300), 500(200), 1000(100), 2000(100), 3000(50), 
5000(10) and 10000(10) from top to bottom, where the inside 
of parentheses represent numbers of realizations for Cp{t; tcq). 
For Cp{t; 0), the number is 1000 for = 1000 and is 100 for 
the others. 



In Stage-I (and III), we fit Cp{t;0) (resp. Cp{t;tcq)) 
by a stretched exponential function, which has three pa- 
rameters: Cp(0;0), tcorr(O) and /3(0) (resp. Cp(0;icq), 
tcoiT(ioq) and /3(toq)). In order to obtain scaling laws for 
the parameters, we show them as functions of degrees of 
freedom N. The parameters Cp(0;0) and Cp(0;teq) rep- 
resent temperature at i = and at equilibrium respec- 
tively, and hence they do not depend on N. We there- 
fore focus on the other four parameters, icorr(O), /3(0), 
icorr(ieq) and /3(ioq)- The correlation functions, Cp(i;0) 
and Cp{t; teq), are shown in Fig llll and values of the four 
parameters are reported as functions of N in Figll2l 

For large N, N > 200, the correlation times are pro- 
portional to N, that is icoir(O) = and icoiT(icq) = 
N/2, and the stretching exponents /?(0) and /3(<eq) are 
almost constants. We expect that these scaling lows for 
the four quantities are kept even in the thermodynamic 
limit, although they break for small N where icorr(O) is 
larger than N/5 and /3(0) is smaller than the constant. 
The duration of Stage-I, ti/u, is around 23 for N = 100, 
and hence icoir(O) and /3(0) are estimated mainly not 
in Stage-I, but in Stage-II from Fig lllT a). In Stage-II, 



icorr(T) and /3(t) are increasing and decreasing functions 
of r respectively (see FigEJ, and hence tcorr(O) and (3{0) 
are larger and smaller than expected values, respectively. 



V. SUMMARY 

As a summary, we have investigated the relation be- 
tween relaxation and diffusion in a Hamiltonian sys- 
tem with long-range interactions. The relaxation pro- 
cess is divided into three stages: quasi-stationary, re- 
laxational and equilibrium. We showed that diffusion 
becomes anomalous only in the second non-stationary 
stage, where magnetization is increasing and goes to- 
wards to the canonical value. The result mentioned above 
does not depend on the number of degrees of freedom, at 
least from N = 100 to 10000. 

The interval where the anomaly in diffusion appears 
becomes shorter and shorter in logarithmic time scale as 
N increase corresponding to a sharper change of mag- 
netization. Moreover, a detailed investigation exhibits 
the absence of power-type diffusion even in the non- 



8 



0.4 
0.3 



n 



Q 



N/2 

.0" 



□ A'/5 



□ 



□ □ 



□ □□ □ □ 



□ 



□ 



□ 



(b)- 



□ □ 



N 



FIG. 12: Two parameters of correlation function as functions 
of degrees of freedom, tcorr{N) (a) and P{N) (b). In both (a) 
and (b), squares (□) represent values of the parameters for 
Cp{t;0) (Stage-I), and crosses (x) for Cp{t;tcq) (Stage-Ill). 



Among the three parameters, Cp(0; r), tcoii{T) and /3(r), 
the stretching exponent /3(r) plays a crucial role to yield 
anomaly in diffusion. If we assume that /3(t) is a con- 
stant, we never observe anomaly in diffusion. This result 
is consistent with the fact that anomaly in diffusion does 
not appear in (quasi-)stationary state, because correla- 
tion function Cp(t;T), and (3{t) accordingly, are almost 
invariant with respect to r. 

We also investigated scaling laws concerning degrees of 
freedom N. The duration of quasi-stationary stage is pro- 
portional to N^'"^ , and relaxation time, at which the sys- 
tem reaches at equilibrium, is also proportional to 7V^'^ 
asymptotically although some corrections must be added. 
In both quasi-stationary and equilibrium stages, ^corr is 
proportional to N, and /3 is almost constant. These sim- 
ple scaling laws imply that fitting by stretched exponen- 
tial functions is valid irrespective of degrees of freedom. 

We have not understood the theoretical reason of 
the appearance of a stretched exponential function. 
If we assume that several time scales with exponen- 
tial correlation function, exp(~t /tcon) , are present, and 
we assume probability distribution function of icorr, 
-P(icorr), then we obtain a stretched exponential function 
/ Pitcorr) expi[— t/tcorr)dtcorr by choosiug Suitable forms 
for P(tcorr) |la|23- In our model, P(tcorr) corresponds 
to the distribution of time scales of individual rotators. 
The investigation of the macrovariable C(t; r) in relation 
with the microvariables of the individual particle corre- 
lation functions will be a subject of future work. 



stationary stage. We guess that anomaly in diffusion is 
a transient anomaly due to non-stationarity. 

Diffusion is obtained by integrating the correlation 
function of momenta Cp(t;T) and the correlation func- 
tion is approximated by a series of stretched exponen- 
tial functions Cp{t;T) = Cp(0;T) exp[-(t/tcorr(T))'^(^)]. 
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